Casimir effect in a one-dimensional gas of free fermions 
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We compute an analog Casimir effect in a one-dimensional spinless Luttinger liquid confined 
to a segment in the presence of a nearly-impenetrable partition dividing the segment into two 
compartments. The Casimir interaction is found to be a bounded piecewise-continuous oscillatory 
function whose maxima are points of force discontinuity and correspond to resonant tunneling across 
the partition. The well-known regularization-based results are reproduced by the lower envelope of 
this function, which corresponds to an approximation that ignores the rather large oscillations due to 
particle discreteness. These macroscopic conclusions are tested and confirmed via a rigorous analysis 
of the Casimir effect in an exactly-solvable model of a one-dimensional non-relativistic spinless gas 
of free fermions, thus resolving an objection that has been raised by Volovik (2003). Additionally 
we confirm the result of a recent calculation which employed an effective low-energy theory with 
a cutoff to find the Casimir interaction between two strong well-separated impurities placed in a 
Luttinger liquid. 

PACS numbers: 71.10.Pm, 73.21.Hb, 03.75.Ss, ll.10.-z 
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I. INTRODUCTION 

Casimir interactions are macroscopic manifestations of 
the vacuum zero-point energy brought about by introduc- 
tion of boundaries or changes in the topology of space. 
These modify the spectrum of zero-point fluctuations 
from their form for free space, which gives rise to exper- 
imentally detectable forces Although originally de- 
rived as an attractive interaction between perfectly con- 
ductive parallel plates induced by vacuum fluctuations 
of the electromagnetic field [2| , similar effects are present 
for other fields, boundary conditions, and geometries [|[. 
Analogous phenomena take place in classical physics due 
to equilibrium thermal fluctuations [J]. 

Although important in their own right, the existence 
of Casimir interactions has a series of implications span- 
ning various areas of science 0. For instance, the van 
der Waals counterparts of these forces are crucial in un- 
derstanding the phenomena of wetting Q . In nanoscale 
devices Casimir attractions may lead to an undesirable 
effect of " stiction" of closely positioned parts of an appa- 
ratus [3] ■ Casimir forces can also be made useful; the first 
micromachines employing them were recently built Q. 
Experimental measurements of the Casimir force provide 
constraints on the existence of extra space dimensions 
and fundamental physics beyond the standard model Q , 
while dynamic versions of the phenomenon [i[ are related 
to the Hawking radiation from black holes and Unruh ef- 
fect of radiation from accelerating objects. 

The Casimir interaction Eq can be defined as a dif- 
ference between the vacuum energy E vac of the system 
constrained by the boundaries or changes in topology and 
the vacuum energy of free space: 
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the case of the electromagnetic Casimir effect), thus im- 
plying that E vac is the sum of zero-point energies of a 
collection of simple harmonic oscillators with a spectrum 
w(k) = c|k|: 
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The vacuum energies, in turn, are calculated from a har- 
monic field theory (such as quantum electrodynamics in 



where the summation is performed over all branches of 
the spectrum and over all allowed wave vectors k. The 
Casimir interaction arises due to a difference between the 
allowed set of wave vectors k of the constrained system 
and that of free space. 

In any realization of the Casimir effect, the boundaries 
are made of real materials that are transparent at ex- 
tremely short wavelengths, so that the high-energy parts 
of Eq. @ are insensitive to the presence of the bound- 
aries thus implying a cancellation in Eq.JT]). For prac- 
tical calculations, one studies the case of an impenetra- 
ble partition and applies a regularization scheme which 
amounts to the introduction of a function /(A, k) such 
that /(A — > 0,k) = w(k). The role of the parameter A 
consists in suppressing the large-k contributions into the 
vacuum energies to guarantee their convergence. Then 
the Casimir interaction is found from Eqs. |T]) and 
by taking the limit A — ► 0. 

For suitable choices of the function /(A, k), the out- 
come of this procedure gives a finite result that depends 
only on Planck's constant, the speed of light and macro- 
scopic length scale(s) indicating that we are dealing with 
a long-wavelength phenomenon; the same result has been 
obtained with different regularization schemes, so that it 
has a universal character. There is general agreement 
on the correctness of the result, and yet the route to it 
is subtle; it is certainly possible to get other (physically 
meaningless) evaluations of the sum @. It is a little 
surprising that in the cancellations that provide a finite 
result for Eq. |T]), no trace of the possible differences 
between long-wavelength and short-wavelength physics 
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survives, even though this would significantly affect the 
sum ©. 

To shed some light on this and other issues having 
their roots in diverging vacuum energy, Volovik [l(j | ad- 
vocated an approach consisting in studying analog ef- 
fects in condensed matter systems. Here the role of the 
physical vacuum is played by the ground state of a con- 
densed matter system, and Casimir-type effects occur 
whenever the ground state is disturbed by boundaries 
or changes of space topology. The ground-state energy 
density of a condensed matter system is finite, and nei- 
ther the harmonic approximation nor regularization are 
needed to calculate the Casimir effect. This is because 
the underlying microscopic physics is governed by the 
well-established quantum mechanics of non-relativistic 
particles interacting via non-singular interactions. As a 
result any quantity of interest can be calculated, at least 
in principle. When the low-energy physics of a candidate 
condensed matter system is described by a harmonic field 
theory, we ought to find the same interaction as is pre- 
dicted by a regularization scheme. 

There are numerous condensed matter systems whose 
low-energy physics is described by harmonic field theo- 
ries. Superfluid 4 He in its ground state is the simplest 
analog to the vacuum of quantum electrodynamics, where 
the role of photons is played by phonons propagating with 
sound velocity c. However the dispersion law u> = c|k| is 
only applicable as long as the phonon wavelength 27r/|k| 
is significantly larger than interatomic distance. 

Volovik provided an example of an exactly-solvable 
toy condensed matter system [ll| where the Casimir 
interaction seemed to disagree with the regularization- 
based conclusion (which is a one-dimensional counterpart 
of Casimir 's original force between perfectly-conducting 
plates). He considered a one-dimensional gas of non- 
relativistic free spinless fermions in the geometry of 
Fig. la, where N particles are placed inside an impen- 
etrable box of length L with a nearly-impenetrable par- 
tition (shown in gray shade). Although this is a non- 
interacting system, nontrivial correlations are built in via 
the Pauli principle. The partition breaks the system into 
two nearly-impenetrable boxes of length a and L — a; an 
infinitesimally small tunnelling transparency allows par- 
ticle exchange between the boxes to minimize the energy. 
The number of fermions M captured between one of the 
walls of the box and the partition can be found by min- 
imizing the ground-state energy E(M, a,N — M,L — a) 
with respect to M under the constraint that M can only 
take on integer values. As a result the M(a) dependence 
has the form of a staircase with plateaus located at inte- 
ger values of M. The inter-plateau transitions are con- 
nected via resonant tunneling since the ground state here 
is degenerate: E(M, a,N — M,L — a) = E(M ±l,a,N — 
M ± 1, L — a). Volovik explicitly demonstrated that at 
the points of degeneracy the force exerted on the parti- 
tion, —dEc /da, undergoes a discontinuity whose mag- 
nitude is parametrically larger than that of the ordinary 
Casimir force implied by the regularization-based result - 
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Figure 1: Geometries employed to study Casimir effect in one 
dimension: (a) box of length L with N interacting fermions or 
repulsive bosons and a nearly-impenetrable partition (shaded 
gray) placed a distance a away from one of the box walls. For 
impenetrable box walls the ground-state energy of this sys- 
tem is identical to that shown in (b): N particles in a ring 
of circumference L with two nearly-impenetrable partitions. 
Regular and bold typefaces show correspondence between ge- 
ometries (a) and (b). 



even in the limit of macroscopic occupancy of both boxes. 
He called this phenomenon the " mesoscopic" Casimir ef- 
fect because "high-energy" physics (arising from the dis- 
creteness of the particles constituting the " vacuum" ) ex- 
plicitly enters into the result even when the sizes of the 
boxes a and L — a are much larger than average inter- 
particle spacing. A similar conclusion was reached for a 
one-dimensional model of " ultra-relativistic" spinless free 
fermions with a linear dispersion law [l2j ■ 

More recently the Casimir interaction between two 
well-separated impurities immersed in a one-dimensional 
spinless quantum liquids was studied within an effective 
low-energy theory with a cutoff [l3[ and in a lattice model 
[l4l |. A closely related problem is that of the interaction 
between two vacancies in the one-dimensional spin-1/2 
antiferromagnetic Heisenberg model These studies 
are related to the L 3> a limit of the ring geometry shown 
in Fig. lb in the presence of two barriers. It is easy to re- 
alize that when the barriers become barely penetrable 
the ground-state energy of any quantum liquid in the ge- 
ometry of Fig. la approaches that in Fig. lb. Thus for 
non-interacting fermions and strong impurities it should 
be possible to compare the conclusions of Refs. [H[ and 
[ilj with Volovik's findings [HI ■ 
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It was found [l3| - [l5| that the Casimir interaction is a 
variable-sign oscillatory function of the impurity /vacancy 
separation a, with maxima for the separations which per- 
mit resonant tunneling. The maxima of the Casimir en- 
ergy represent the points of discontinuity of the slope. 
Qualitatively this agrees with Volovik's results [llj]. We 
also note that an oscillatory interaction between vari- 
ous kinds of inhomogeneities immersed in an ideal three- 
dimensional Fermi gas was found by Bulgac et al. 

The goal of this paper is two-fold: first contributions of 
the effects of the particle discreteness into Casimir-type 
phenomena are computed for an arbitrary spinless quan- 
tum liquid by employing the conventional macroscopic 
field-theoretical viewpoint. Second, these phenomeno- 
logical conclusions are verified through a rigorous mi- 
croscopic analysis of the Casimir effect for spinless free 
fermions. As a result we will be able to: (i) establish a re- 
lationship between the ordinary and mesoscopic Casimir 
effects, and (ii) to link the macroscopic and microscopic 
approaches. 

In addition to these issues of principle, the prob- 
lem considered here may be of practical interest as the 
Casimir forces between impurities immersed in fermionic 
quantum wires realized with ultracold gases [l7| might 
be experimentally measured in the future. 

Before proceeding we hasten to mention that the ef- 
fects calculated in this work should not be confused with 
the relativistic fermionic Casimir effect (with bag-model 
boundary conditions on the surface), which for parallel 
plates was first worked out by Johnson [l8j . 

The organization of this paper is as follows. In Sec- 
tion II we summarize conventional predictions of the har- 
monic theory and explain how the complications com- 
ing from particle discreteness can be accounted for phe- 
nomenologically within the harmonic theory. Section III 
turns to the non-interacting fermion system, showing how 
these results can be verified and confirmed via rigorous 
microscopic analysis. Section III A deals with the case of 
a segment whose ends are subject to either an arbitrary 
combination of the Dirichlet and Neumann or to periodic 
boundary conditions imposed on the many-body wave- 
function. Here the Casimir effect manifests itself as a 
finite-size correction to the ground-state energy. The ir- 
relevance of the effects of the particle discreteness in the 
segment or ring geometries allows us to probe only the 
continuum part of the Casimir force. 

A more complicated geometry where the effects 
of the particle discreteness come into play (a seg- 
ment with nearly-impenetrable partition or two nearly- 
impenetrable partitions inserted into a ring) is studied in 
Sections IIIB and IIIC. The former deals with the limit 
when the particle discreteness continues to be ignored. 
In close correspondence with the conclusions of Section 
IIIA, we reproduce the results of the effective theory for 
the case of the Dirichlet boundary conditions enforced at 
the segment ends (which is equivalent to the case of two 
partitions inside a ring). The effects of particle discrete- 
ness are incorporated in Section IIIC where macroscopic 



results of Section IIA are confirmed. We conclude (Sec- 
tion IV) by discussing the correspondence between mi- 
croscopic and macroscopic viewpoints and posing several 
questions for future study. 



II. MACROSCOPIC APPROACH 

Volovik's fermion model exhibits large oscillations in 
the Casimir energy, which arise from effects of particle 
discreteness. In the following subsections we show that 
these can be understood phenomenologically without re- 
sorting to microscopic analysis. 

A. Continuum results 

We consider a generic harmonic spinless quantum liq- 
uid in the geometry of Fig.l. We begin with restat- 
ing the well-known result [19| that the long-wavelength 
low-energy dynamics of a translationally-invariant one- 
dimensional spinless quantum liquid comprised of inter- 
acting fermions or repulsive bosons of mass m is governed 
by the harmonic (Luttinger liquid) theory with the La- 
grangian 



( du ]2 
K dt' 
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where x is the position, t is the time, and n is the num- 
ber density. The field u(x,t) represents displacement 
of the particles with respect to the case of equidistant 
positions. The first term of describes the kinetic 
energy while the second accumulates pertinent informa- 
tion about the interparticle interactions (and statistics, 
in case of fermions) parameterized by the sound veloc- 
ity c. This latter quantity is determined by the ground- 
state energy; for spinless free fermions the sound velocity 
c = Trhn/m is identical to the Fermi velocity. The de- 
scription §5§ is applicable for spatial scales exceeding the 
interparticle spacing n , provided that the underlying 
interactions are sufficiently short-ranged and the strain 
is small (\du/dx\ <§; 1). 

The interactions due to finite size described above will 
have the universal Casimir form but the boundary condi- 
tions also matter. In cases of interest the collective vari- 
able u(x,t) satisfies either the Dirichlet condition u = 0, 
or the Neumann condition du/dx = 0. The former mim- 
ics the perfectly-conducting plate of the original Casimir 
calculation while the latter corresponds to a plate of 
infinite magnetic permeability [20| . If in the geometry of 
Fig. la the partition and the walls enforce "like" bound- 
ary conditions (i. e. either Dirichlet or Neumann at all 
boundaries), then for any smooth regularization scheme 
(and up to a constant) the Casimir interaction was found 
to be of the form [201 
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The case of mixed boundary conditions is particularly 
interesting because Boyer demonstrated |2(| that f° r the 
geometry of Fig. la the Casimir interaction is repulsive 

ttHc f 1 1 \ 

Ec,unlike = h 7 (5) 

4e \ a L — a J 

when the partition and the walls enforce " unlike" bound- 
ary conditions 2(j (i.e. either Neumann at the walls and 
Dirichlet at the partition or Dirichlet at the walls and 
Neumann at the partition). 

For L — > 00 Eq.jH gives Ec.uke = — 7rfic/(24a) which 
is a well-known result [H[, 0, H3, while E q-© pre- 
dicts Ec\ U niike = irhc/ (48a). Regardless of the boundary 
conditions employed, the he/ a behavior of the Casimir 
interaction Ec can be understood via dimensional analy- 
sis. Indeed, since the effect is both quantum-mechanical 
and " relativistic" , Ec must depend on a, h and c. If 
we additionally assume that the effect is universal, then 
Ec cannot depend on microscopic quantities such as the 
mass of the underlying particles or the wavevector cut- 
off n. Given Ec, h, c and a one can form one and 
only one dimensionless combination Eca/hc which im- 
plies Ec = const (He /a) where the numerical constant is 
only determined by the boundary conditions. 

B. Effects of particle discreteness 

The Lagrangian ^ describes the elastic theory of a 
one-dimensional solid. This viewpoint is acceptable be- 
cause despite the fact that the long-range translational 
order is destroyed by the zero-point motion [23[ , algebraic 
correlations survive and typical strains remain small. 
Quantum fluctuations are additionally suppressed under 
confinement as in geometries of Fig.l. 

A complication arises in the discussion of the elastic 
theory of a solid, in which the field variable describes 
the displacements of discrete particles from equilibrium 
positions. When we change the length of a box con- 
taining a certain number of particles, we change both 
the wavelengths of the excitations of the system and the 
equilibrium positions themselves: the solid experiences 
a strain, which changes the ground-state energy. When 
the strain becomes large enough, it can be relaxed (as 
in geometries of Fig.l) by adding or removing a particle, 
but this is a discrete effect. A second complication is 
that the introduction of the boundaries will cause parti- 
cle displacements, no matter where they are put, because 
the particles might be attracted or repelled by the bound- 
aries. These effects can be described within the harmonic 
theory ([3]), as we will now demonstrate. 

The ground-state energy of a quantum liquid in the 
geometry of Fig. la can be obtained by minimizing the 
sum of the ground-state energies of separate arbitrarily 
populated boxes E(M, a) + E(N — M, L~a) with respect 
to M provided the latter is restricted to the integers. In 
the macroscopic limit M, N — M 3> 1 when the discrete- 
ness of the underlying particles can be at first neglected, 



this condition reduces to the statement of the equality of 
the chemical potentials on the two sides of the partition. 
This in turn implies that the bulk densities (and thus 
the sound velocities) on the two sides of the partition 
are the same. The common density n defines the refer- 
ence zero-strain state, and the field variable u in Eq.((3|) 
describes displacements away from corresponding inter- 
particle spacing n~ x . Although the chemical potentials 
on the two sides of the partition are equal, the forces are 
not since the partition is not in mechanical equilibrium. 
The corresponding osmotic force on the partition is the 
continuum Casimir force implied by Eq.([J| and ([5]). 

When the discreteness of the particles is taken into 
consideration, a particle exchange across the partition 
can only occur at the points of degeneracy of the ground 
state. As explained above generally this causes a strain 
which can be computed as follows. The number of par- 
ticles captured between x — and x — a, namely 
n(a — u(a,t) + u(0,t)) — (p, is constrained to be an in- 
teger M. This leads to the boundary condition 

u(a, t) — u(0, t) = a (6) 

n 

Here <p is a correction of order unity due to the distur- 
bance of the density profile away from uniformity caused 
by the ends of the interval. This correction cannot be 
computed within the harmonic approach since the latter 
fails in the vicinity of the left wall and the partition. 

The boundary condition ([5]) implies that in general 
there will be a uniform static strain of the form 

du{x)/dx = "(M)-"(0,*) = 1 M + v 
a na 

The system can also have time-dependent motions in ad- 
dition to the static strain; these are described by the wave 
equation d 2 u/dt 2 — c 2 d 2 u/dx 2 = 0. 

For example, if the field u(x, t) obeys the Dirichlet 
boundary condition at the x — wall, u = 0, then for 
^ x ^ a the corresponding general solution to the wave 
equation satisfying the boundary condition ^ is given 
by 

u ( Xj t ) = (1 - ^ + *)x + Y bl e Mct / a sin — (8) 
na A — ' a 

i=i 

where the b\ are arbitrary complex coefficients. The sum 
in Eq.((5|) represents the field configurations responsible 
for the continuum Casimir effect in an interval of length 
a with Dirichlet-Dirichlet (DD) boundary conditions im- 
posed at the segment ends. Substituting the solution ([8]) 
into the Lagrangian ([3]) we find 

C = tf DD )-%£(na-M-<p) 2 (9) 
2na 

where C^ DD ^ is the Lagrangian corresponding to the sum 
in Eq.©. It is important to notice that the cross terms 
occurring in the evaluation of (du/dx) 2 vanish upon inte- 
gration over x. This implies that the term linear in x and 
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the sum in ([8]) contribute independently to the ground- 
state energy. Therefore the portion of the Casimir energy 
coming from the ^ x ^ a modes is given by 



E (DD) 



+ - — (na-M-wY 
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irhc 
Ma' 
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(na-M-ip) 2 
(10) 

where E^ DD > = -Tthc/(2Aa) is the first term of Eq.g]). 
The dimensionless parameter 



mc 



(11) 



in the second term of (|10[) can be recognized as the Lut- 
tinger liquid exponent governing the low-energy behavior 
of the correlation functions in one-dimensional harmonic 
liquids 19]; the free-fermion case corresponds to g = 1. 

For the geometry of Fig. la there are two segments. 
The portion of the Casimir energy coming from the 
modes in the segment a < x ^ L is computed similarly. 
Combining it with Eq. (|10[) we arrive at the expression 



E, 
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where we specified to the case of the Dirichlet wall at 
x = L (which explains the third "D" in the "DDD" su- 
perscript) and employed the continuity of the field vari- 
able across the partition at x = a. Since the latter is 
slightly transparent to particles, the integer parameter 
M can adjust to minimize this expression for any chosen 
value of a. The result is that the ground-state energy of 
the system is an oscillatory function of the position of the 
partition. The usual Casimir effect is the lower envelope 
of the function, while particle discreteness contributes 
parabolic cusps. It should be noted that the scale of the 
discreteness effect is comparable to the contribution from 
the zero-point motions. The conclusion (|12p is also ap- 
plicable to the geometry of Fig. lb; in the a <C L — a limit 
it reduces to the Casimir interaction between two strong 
impurities placed in an infinite Luttinger liquid [T3 |. 



III. MICROSCOPIC ANALYSIS OF SPINLESS 
FREE FERMIONS 

In the following subsections we will calculate the 
ground-state energy for a model of non-interacting spin- 
less fermions. Since the long-wavelength excitations of 
this system are described by a harmonic theory, the 
finite-size (Casmir) effects are captured by the contin- 
uum theory as described in previous Section. Here we 
will be able to see how these macroscopic conclusions 
emerge from a completely different microscopic analyis of 
the exactly-solvable model. Naturally this will not have 
the generality of phenomenological theory, however all 
the details, including the phase of the Casimir oscillation 
ip in Eq. (112|) will be determined. 



There is again an issue of boundary conditions imposed 
on the many-body wavefunction (hereafter indicated by 
the letters D or N) . The (D) Dirichlet boundary con- 
dition for the fermion wavefunction implies that no par- 
ticles can reach the boundary which is consistent with 
the Dirichlet condition imposed on the collective variable 
u(x,t) (i.e. u — at the boundary). However, the Neu- 
mann condition for the collective variable (du/dx = 0) 
allows a current through the boundary, while the (N) 
Neumann condition for the wavefunction does not. At 
the same time the (N) Neumann condition for the wave- 
function is not inconsistent with the Dirichlet condition 
imposed on the collective variable. These statements rep- 
resenting plausible guesses will be verified below; they im- 
ply that the conclusions derived from microscopic analy- 
sis of the Dirichlet and Neumann geometries can be only 
compared with the macroscopic Dirichlet results (|4|) and 

G3- 



A. Segment and ring geometries 

Before directly addressing the geometries of Fig.l we 
find it useful to start with simpler segment and ring ar- 
rangements which may be viewed as reference states for 
those in Fig. 1. Here the Casimir interaction can be ex- 
tracted as a finite-size correction to the bulk ground-state 
energy [221 ]. 

It is intuitively clear and will be made explicit below 
that the rearrangements of the density profile caused by 
either Dirichlet or Neumann boundary condition repre- 
sent strong perturbations which cannot be correctly cap- 
tured within the harmonic theory ([3]). 



1. Like boundary conditions 

First consider M fermions inside an impenetrable box 
of length a. This corresponds to the "D" boundary con- 
ditions enforced at the segment ends. The normalized 
single-particle wavefunctions for this problem are well- 
known: 



1/2 



sin kix 



(13) 



where the allowed wave numbers satisfy the condition 
ki = nl/a, I = 1,2,3,.... The ground state is built by 
sequential occupation of the single-particle states begin- 
ning from I — 1 and ending with I — M . The resulting 
expressions for the ground-state energy 
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and the particle density inside the box 



M 



Pdd{x) =^2ipf(x) 



M 



1=1 



sin (2ir(M + \)x/a) 
2a sm(Trx/a) 

(15) 

are well-known [24j . Since the latter is depleted near 
the box walls, the average bulk density, as evidenced by 
Eq. (fTS"|) . is (M + l/2)/a rather than the naively expected 
M/a. Swiatecki [251 ] emphasized the importance of this 
fact in separating the bulk contribution to the ground- 
state energy from those due to the boundaries. We made 
this explicit in Eq. (fl"4")) which has a form typical of the 
ground-state energy of a one-dimensional quantum sys- 
tem confined to a segment of length a 22]: The first 
term gives the bulk energy of the spinless free Fermi 
gas of density n = [M + 1/2) /a while the second term, 
—Tr 2 h 2 n/(2Ama) = — 7rfi,c/(24a), represents the Casimir 
effect 

If both ends of the segment enforce "N" boundary 
conditions, the single-particle ground-state has zero en- 
ergy and is characterized by a constant wave function, 
ipo{x) = l/y/a. The excited states are of the form, 
ipl(x) — (2/a) 1 / 2 coskix, and the allowed wave numbers 
are determined by k\ — nl/a, I — 1,2,3,.... It is then 
straightforward to realize that the ground-state energy 
of the M-fermion segment follows from Eq. (Ti"4")) by sub- 
stituting M — 1 for M: 
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The particle density is then given by 



Pnn(x) 



M-l 
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1=0 



M-\ sin(27rx(M- |)/o) 
a 2asm(Trx/a) 

(17) 

Eq. p?|) makes obvious that fact that " N" boundary con- 
ditions promote particle accumulation near the segment 
ends where the density reaches the maximal value of 
(2M — I) /a. This is to be expected as the "N" boundary 
condition can be implemented via an attractive potential 
localized at the boundary. As a result the bulk density 
n = (M — 1/2) I a is smaller than M/a. This is con- 
sistent with the expression for the ground-state energy 
(|16p which implies an attractive Casimir interaction of 
the expected — 7r?ic/(24a) form [221 ]. 

Although Eas.([i"4]) and (fl"6|) do not explicitly address 
the geometry of Fig. la to which Eq.lfJJ is supposed to 
be applicable, they certainly imply plausibility of Eq.(fi]) 
for like boundary conditions. 



single-particle wave functions are still given by Eq. (|13|) 
while the wave numbers selected by the "N" boundary 
condition at x — a satisfy fc; = ir(l— l/2)/a, I = 1,2, 3, .... 
As a result the ground-state energy and the particle den- 
sity inside the segment are given by 
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and 



Pdn{x) 
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(19) 



respectively. The latter vanishes at the "D" end x = 
and reaches the maximal value of 2M/a at the "N" end 
x = a. As evidenced by Eq. (fT9)) . the particles pushed 
away from the "D" end accumulate at the "N" end, and 
as a result the bulk density n takes on the naively ex- 
pected M/a value. Then Eq. (fT8"]) implies an attractive 
Casimir interaction, —irhc/(24a), of exactly the same 
form as was found for like boundary conditions (see 
Eqs.(H} and HID). 



3. Periodic boundary conditions 

Now assume that M fermions are confined to a circle 
of circumference a. For odd M the single-particle wave- 
functions for this problem are of the expifcjx form with 
the wavenumbers satisfying the condition ki — 27rZ/a, I — 
Q,±l,±2, .... The ground state is built by placing one 
of the particles into the I = state while the remaining 
M — 1 fermions symmetrically occupy the states begin- 
ning from I = ±1 and ending with I — ±(M — l)/2. The 
ground-state energy is thus given by 



Eperiodic(-M , a) 



n 2 h 2 (M_ 
6m V a 



ir 2 h 2 M 
6ma a 



(20) 



Repeating the construction for even M leads to the same 
result. Eq. ([20]) has a structure similar to Eqs. tfT!"]) , 
(fTB)) and (|18p : the first term gives the bulk energy 
of spinless gas of free fermions of density n — M/a 
(which is everywhere uniform) while the second term, 
—ir 2 h 2 n/(6ma) = —irhc/ (<oa), represents the Casimir ef- 
fect. The magnitude of the former in the ring geometry 
is four times larger than that for the segment [22[ ■ 



B. Geometries with three boundaries; continuum 
limit 



2. DN boundary conditions 

As a next step consider M fermions belonging to a 
segment of length a with the " D" boundary condition at 
x = and "N" boundary condition at x = a. Then the 



The geometries of Fig. 1 may be viewed as a result 
of insertion of one (a) or two (b) partitions inside a seg- 
ment or a ring of length L containing N fermions. This 
brings up a qualitatively new aspect into the problem: 
the issue of particle discreteness. This was not relevant 
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to the segment or ring geometries since the total particle 
number was fixed. Now as the partition (or partitions) is 
adiabatically displaced, particle exchange across it has to 
be allowed - otherwise instead of the Casimir effect one 
would be probing the bulk compressibility of the Fermi 
gas. The effects of particle discreteness are largest when 
the partition is nearly-impenetrable, i. e. it enforces the 
"D" boundary condition. In what follows we will only 
focus on this case. However for the geometry of Fig. la 
the boundary conditions at the segment ends could be 
either like (DD or NN) or unlike (DN or ND). 

We begin with the case that the partition(s) and the 
walls enforce Dirichlet constraints at all three boundaries. 
When the number of particles captured in each of the 
boxes in Fig.l is large, these may be treated as continuous 
varibles. The equilibrium value of a macroscopic variable 
(for example M) can be found by minimizing the ground- 
state energy, i. e. by requiring that d(Ez)z)(M, a) + 
E DD (N — M,L — a))/dM = 0. It is straightforward to 
verify that to leading order in M , N — M 3> 1 this gives 



M+i _N-M+ \ _N- 



1 



(21) 



which, in view of Eq. (|15|> . can be recognized as the 
condition that the bulk densities in the two boxes be 
equal. Their common value of (N + l)/L is higher than 
(N+ 1 /2) jL of the reference iV-fermion segment of length 
L or N / L for the ./V-fermion ring of circumference L be- 
cause the partition(s) push the particles into the bulk. 

Combining Eqs. (jT4j) and ([21] ) we find that the 
total energy for the geometries of Fig. 1, Eddd = 
E DD (M, a) + E DD {N - M, L — a), is given by 



2 ^2,^,3 



E 



k 2 K 



DDD 



6 m 



-L 



nhc ( 1 



1 



24 V a L 



(22) 



where n = (N +1)/L and we used the fact that irh(N + 
l)/(mL) = irhn/m = c is the sound velocity. 

We observe that the first term of Eq. (f2"2"|) is of or- 
der O(L) and represents the bulk energy of a Fermi 
gas of density n = (N + l)/L. The remaining in- 
dependent terms combined which are of order 0(a _1 ) 
and 0((L — a) -1 ) are identical to the result based on 
regularization, thus supporting Eq.(U]). Since the O(l) 
term is absent, we conclude that there is zero self-energy 
associated with a single partition. This is consistent with 
Eq. (|14[) which also predicts zero self-energy but in vari- 
ance with the earlier conclusions of Swiatecki [25| and 
Peierls [Hj]. 

The analysis of the NDN and DDN geometries (Fig. 
la) leads to the same result (|2"2"j) with n — N/L and 
n = (N + \)/L, respectively. 



C. Geometries with three boundaries; effects of 
particle discreteness 

The main conclusion of the previous analysis is that 
when the particle discreteness is ignored, the Casimir 



interaction is universally given by Eq.©. Specifically, 
for the segment geometry of Fig. la it is insensitive to 
the boundary conditions enforced at the segment ends 
(provided they form any combination of "N" and/or "D" 
constraints). As a next step we improve the analysis by 
taking into consideration the particle discreteness. 



1. DDD boundary conditions: nearly-impenetrable 
partition(s) and "D" walls, Fig. 1 

In order to capture only the a-dependent contributions 
the Casimir energy for this case will be defined as 

E c ,ddd = E DD (M,a) + E DD (N-M,L-a) 
n 2 h 2 (N + lf 



6mL 2 



(23) 



and it is assumed that for given a an integer M is cho- 
sen to minimize (f23| . We already demonstrated (see 
Eq. (f2"2"|) ) that if the discrete character of M is ignored 
(and M,N — M > 1) then E c ,ddd reduces to Eq.®. 
Since the latter is the result of unconstrained minimiza- 
tion with respect to M, while M must be an integer in 
the fermion model, we conclude that the effect cannot 
be smaller than its continuum approximation. In other 
words, the continuum approximation is the lower enve- 
lope of the interaction given by Eq. (|2"3"]) . 

Microscopic consideration. In order to understand 
the Ec,ddd(o) dependence (shown in Fig. 2) , let us begin 
with a very small size a for one of the boxes in Fig. 1. 
Then the energy cost of localizing even one particle inside 
the box is very large; as a result this box is empty and 
all the particles reside inside the second box of size L — a. 
This continues to be the case as a adiabatically increases 
as long as the latter remains sufficiently small. Thus the 
Ec,ddd(o>) dependence is given by the M = branch of 
Eq. (f23|) - it is an increasing function of a reflecting the 
compression of N fermions captured inside the box of 
decreasing size L — a. As a — > the force exerted on the 
partition (Fig. la) or the interaction force between two 
partitions (Fig. lb), F = —dEc,DDD/9a, approaches 
-Tr 2 h 2 (2N 3 + 3N 2 + N)/(6mL 3 ) which remains finite in 
the limit that N, L -> oo with N/L fixed. The M = 
regime obtains until the empty box becomes so large that 
the energy cost of localizing a fermion inside it equals 
the strain energy released when a particle tunnels out 
of the compressed box. At that instant, Eud(N,L — 
a) = E DD (l,a) + E DD (N -1,L- a), the E c , DDD {a) 
dependence switches to the M = 1 branch of Eq. (|23| . 

As a continues to increase, the function Ec,ddd(<i) 
first decreases because the associated falloff of the energy 
of localization of a fermion inside the box of size a prevails 
over the increase of energy of the compressed (N — 1)- 
fermion box. This trend however is eventually reversed 
and the energy cost of compression begins to dominate 
over the energy gain of derealization. As a result, the 
Ec,ddd{o) dependence has a minimum. 
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Upon further increase of a, a point is reached when the 
single-fermion box is too large, while its (N — l)-fermion 
counterpart is too compressed, so that it is energetically 
favorable for a particle to tunnel into the box of size 
a. At that moment, Edd(1, a) + Edd(N — 1, L — a) = 
E DD (2, a) + E DD (N -2,L- a), the Ec,ddd(o) depen- 
dence switches to the M = 2 branch of Eq. (|2"5|) . The com- 
petition between derealization and compression period- 
ically resolved by tunneling keeps recurring as a further 
increases, and the function Ec,ddd(o-) switches through 
the M = 3,4, ...,7V branches of Eq.ipg]), At the transi- 
tion points between the neighboring branches the ground 
state is degenerate, Edd(M, a) + Edd(N — M, L — a) = 
E DD (M + l,a) + E DD (N — M — 1,L — a). All but 
M = and M = N branches of Ec,ddd{o) have min- 
ima where the Casimir force vanishes; their loci should 
be well-approximated by the continuum result ([4]). 

To summarize, the Casimir interaction (|23p is an oscil- 
latory piecewise-continuous function of a whose maxima, 
the points of force discontinuity, correspond to resonant 
tunneling between the boxes in Fig . 1. This agrees with 
previous conclusions [ll|, [l3|, 0, EH • 

Macroscopic limit. When there are large num- 
bers of particles on either side of the partition, the 
function Ec,ddd{«) can be accurately characterized by 
mathematically simpler expressions. The minima of the 
Casimir energy Eq. (|23p occur when the densities on the 
two sides of the partition are equal. This is the condition 
already given in Eq. (|2ip , which can be restated to imply 
that the loci of the zeros of the Casimir force, measured in 
units of the bulk interparticle spacing n^ 1 — Lj (JV + 1), 
are half- integers: 



0.2 



N+l 



cim = naM = M + 



(24) 



where M is the integer that labels the M-th minimum of 
the Casimir interaction. 

The transition points, cim^m+i, between the M-th 
and (M + l)-th branches of the function (j23|) follow 
from the condition of degeneracy, Edd (M, o)+Ed d{N— 
M, L-o) = E DD (M + 1, a) + E DD {N — M — 1,L — a), 
Taylor expanded up to second order in AM — 1. We 
find that in units of the bulk interparticle spacing, the 
transition points are given by 



N + l 
L 



um^m+i = na,M^M+i = M +1, 



(25) 



i. e. they are at integer values. Not surprisingly, the 
distance between the nearest minima of the Casimir in- 
teraction (|2"4"|) is the same as the distance between near- 
est maxima (|25l) and equal the bulk interparticle spacing. 
At the same time the distance between a transition point 
and two nearest minima of (|23p is half interparticle 
spacing. 

Since the loci of the minima of the Casimir interaction 
(|24p are symmetrically sandwiched between the transi- 
tion points (|2l))) . the M-th branch of the function 



Exact 

Continuous 
Macroscopic 




Figure 2: (Color online) Casimir interaction for "D" bound- 
ary conditions enforced at all boundaries, Eq. (|23[) . and its 
continuum, Eq.@ and macroscopic, Eq, (l26[l . approxima- 
tions in units of the bulk Fermi energy e_F = Tr 2 h 2 n 2 /2m, 
as a function of dimensionless position of the partition na 
(Fig. la) or distance between the partitions (Fig. lb) in the 
^ na ^ (N + l)/2 range for iV = 24 fermions. The macro- 
scopic theory is indistinguishable from the exact result, except 
at very small distances. 



may be approximated by its second order Taylor expan- 
sion around a = clm' 



E, 



C.DDD 



(a) 



n 2 h 2 n 2 ( 1 



1 



24m VM+i N-M+\ 



\ ! -1 + \2{na — M - -)' 



nhc ( 1 1 



24 V a L - a 



x ( 1 - 12(na — M - -) 2 



(26) 



which is valid when na is positioned between nearest in- 
tegers 



M < na < M + 1 



(27) 



In the second representation of Eq. ([26| the difference 
between <zm and a was neglected in the coefficients of the 
Taylor expansion which is a valid approximation within 
the range ([27]) . This allows us to see that Eq. ((26|) is 
identical to the macroscopic result (|12p with g = 1 (non- 
interacting fermions) and tp — 1/2. 

As implied by Eg . (|26]) . the oscillations of the Casimir 
force F = — 8Ec : ddd / da are confined between upper 
and lower envelopes whose macroscopic limit is given by 



Eenvelope — i" 



2m 



1 



L — a 



(28) 



The amplitude of this oscillation, 2\F enve i ope \, is the mag- 
nitude of the jump of the Casimir force at the point of 
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resonant tunneling across the partition(s). This result 
was first given by Volovik [Tl| . We should compare this 
with the magnitude of the force computed from the con- 
tinuum theory: according to Eq. it varies as 1/a 2 . 
Then in the 1 <C M <C N limit, the force exerted on 
the partition (Fig. fa) or the interaction force between 
two partitions (Fig. lb) is a factor of M larger than its 
continuum counterpart [11| . 

For the purposes of comparing of the continuum, 
Eq.Q, and macroscopic, Eq. ([2"u) . approximations to 
the exact Casimir interaction, Eq. (|23p . it is natural to 
measure the energy in units of the bulk Fermi energy 
€p = Tr 2 h 2 n 2 /2m, while the position of the partition a 
(Fig. la) or the distance between the partitions (Fig. lb) 
is measured in units of the bulk interparticle spacing n _1 . 
Then the variable na varies between zero and N + I; in 
view of the reflection symmetry about a — L/2 we only 
need to look at the ^ na ^ (N+l)/2 range. The result 
is plotted in Fig. 2 for the case N — 24. 

We note that even for the small-M branches of the 
unapproximated interaction (|23|) the lower envelope is 
well-appproximated by the continuum result ([4|). The 
macroscopic result fEq. (|26p ) represents an excellent ap- 
proximation to the interaction (|23| . Although the M = 
branch is not reproduced, the M = 1 branch is approx- 
imated fairly well, and Eqs. ([23| and Eq. ([26|) are hardly 
distinguishable for M ^ 3. We also note that for an even 
number of particles N there is a minimum of the Casimir 
interaction in the middle of the system at na = (N+l)/2, 
as demonstrated in Fig. 2. Had we chosen N odd, there 
would be instead a maximum at na = (N + l)/2. 

To judge the accuracy of the various approximations 
it is more relevant to look at the experimentally mea- 
surable Casimir force, F = —dEc,DDD / da, rather than 
the interaction. The former (measured in units of nep) 
is plotted as a function of na in Fig. 3 with the same 
choice of N = 24 fcrmions as in Fig. 2. We now see that 
the continuum limit gives a very poor representation of 
the Casimir force. The macroscopic approximation, on 
the other hand, continues to be very accurate; for M ^ 6 
it is hardly possible to distinguish it from the unapprox- 
imated result. The large M (large na) part of the force 
plot resembles Figure 29.3 of Volovik hj. 

As anticipated in Section IIB, the deviation of the 
exact fermionic Casimir interaction from its continuum 
counterpart (J3J) which is due the combined effect of par- 
ticle discreteness and near impenetrability of the parti- 
tion^) can be understood in terms of the effective bound- 
ary condition at the partition(s) ([6|) imposed on the dis- 
placement field u(x, t) of the continuum theory (J3J). For 
the DDD geometry this condition becomes 

u(a,t) = a- M+2 (29) 
n 

where we used ip — 1/2. At the values of a for which un- 
constrained minimization of the ground-state energy with 
respect to M happens to give an integer M, the boundary 
condition in question is "impenetrable". These a's corre- 



1.5 




1.5 3 4.5 6 7.5 9 10.5 12 
na 

Figure 3: (Color online) Casimir force F — — 0Ec,ddd /do, 
and its approximations corresponding to the energy curves of 
Fig. 2. The force is measured in units of neF. Additionally, 
the macroscopic limit of the upper and lower envelopes of 
the force, Eq. (l28|l . is shown. The "macroscopic" and "exact" 
curves almost entirely coincide, and give a force of oscillating 
magnitude and sign that is generally large compared to the 
continuum result. 



spond to the minima of the Casimir interaction where the 
continuum approximation is accurate. Indeed combining 
Eqs. ljM)) and (f2"9")l we arrive at the Dirichlct condition 
u(a,t) = 0. At any other value of a the effective bound- 
ary condition interpolates between the "impenetrable" 
and "transparent" limits. The latter corresponds to half- 
integer nu(a,t) as can be seen by combining Eqs. (|25|) 
and (|29p. The "transparency" is perfect at the points of 
resonant tunneling which are the maxima of the Casimir 
interaction; in the macroscopic limit their loci are im- 
plicit in Eq.(H5]): 



Fenvelope — n I F z ) (30) 

12 \a L — a I 



The latter has the same sign as Eq.© describing the case 
of the Neumann partition (du/dx — 0) and the Dirichlct 
walls (u = 0) but a different magnitude. It is curious that 
the average of the lower (g]) and upper ([3TJ]) envelopes of 
the interaction is given precisely by Eq.©. 

The analysis of the Casimir interaction for the cases of 
"N" or mixed boundary conditions enforced at the seg- 
ment ends is very similar to the one just presented; only a 
summary of the main results emphasizing the differences 
between various cases is provided below. 
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2. NDN boundary conditions: nearly-impenetrable 
partition, "N" walls, Fig. la 

The a-dependent part of the Casimir interaction is 
given by 



ndn = E nd (M, a 



■E DN (N-M,L-a)- 



TT 2 fl 2 N 3 



6mL 2 
(31) 

In the macroscopic limit the loci of the zeros of the 
Casimir force measured in units of the bulk interparti- 
cle spacing -nT 1 = L/N are integer: 



N 



-clm = nau = M 



(32) 



In the same dimensionless units the transition points be- 
tween the M-th and (M + l)-th branches of the function 
(|3ip or, equivalently, the maxima of the Casimir interac- 
tion are given by half- integers: 



N 
T 



«M^M+l = TldM-tM+l 



M 



(33) 



The M-th branch of the function (|3ip may be approx- 
imated by its second order Taylor expansion around 
a = a M : 



E. 



(a) 

C.NDN 



(a) 



n 2 H 2 n 2 (\_ 
~ 24m \~M + 
x (-l + 12(na- 

Tthc ( 1 

~ ~~24~ la + i 



1 



N-M 
M) 2 ) 
1 



x (l - 12(na- M) 2 ) 



(34) 



which is valid when na is positioned between nearest half- 
integers 



M--^na<M + - 
2 2 



(35) 



Eqs. (|3"Tj) - ([331) are counterparts of the set of Eqs. (f23l) - 
(l27f describing the macroscopic limit of the Casimir ef- 
fect with all the boundaries subject to the "D" condition. 
In the macroscopic limit, the only difference between the 
two is the location of the extrema of the Casimir inter- 
action: now the minima lying at integer values of na are 
sandwiched between half-integer maxima of the Casimir 
energy. This is a manifestation of the long-range effect 
exerted by the ends of the segment. We also note that 
Eq. ([34]) coincides with the g = 1, <p = case of the phe- 
nomenological result (fl2"|) . 

These conclusions are illustrated in Fig. 4 where the 
Casimir interaction, Eq. (|3ip . is plotted in dimensionless 
units of the Fermi energy ep as a function of the dimen- 
sionless position of the partition na (Fig. la). Since now 
the bulk density is given by n = N/L, the variable na 
varies between zero and N; in view of the reflection sym- 
metry about a — L/2, only the ^ na ^ N/2 range 




Figure 4: (Color online) Casimir interaction for nearly- 
impenetrable partition and "N" boundary conditions enforced 
at the segment ends, Eq. (|3ip . in units of bulk Fermi energy 
£f as a function of dimensionless position of the partition na 
(Fig. la) in the < na < N/2 range for JV = 24 fermions. 
The inset shows the plot of the corresponding Casimir force. 



is shown. The inset shows the Casimir force. To facili- 
tate comparison with the previously studied case of " D" 
boundary conditions enforced at all boundaries (Figs. 2 
and 3), we again show the case N = 24. It now be- 
comes apparent that the bulk behavior is identical in 
the two cases, apart from a shift by a half along the 
na axis. Additional differences overlooked in the macro- 
scopic approximation are limited to the vicinity of small 
na; they represent manifestations of the distinction be- 
tween the "D" and "N" wall boundary conditions en- 
forced at na = 0. 



3. DDN boundary conditions: nearly-impenetrable 
partition, unlike boundary conditions at the walls, Fig. la 

The a-dependent part of the Casimir interaction is 
given by 



E. 



(a) 

C.DDN 



E D D [M, a)+E DN (N -M,L-a) 



6mL 2 



(36) 



This case is extremely similar to that of all three bound- 
aries constrained to the "D" boundary condition (Section 
IIIC1). In the macroscopic limit it is impossible to tell 
the two apart even when the effects of particle discrete- 
ness are accounted for. The one slight difference (apart 
from the value of the bulk density n = (N + 1/2) / L) 
is that the AT — M + 1/2 combination on the first line 
of Eq.® has to be replaced with N - M. This differ- 
ence however disappears in the second representation of 
Eq. (|2"6")) . The macroscopic limit again coincides with the 
g = 1, (p = 1/2 case of the phenomenological result CC 



11 




Figure 5: (Color online) Top: Casimir interaction for nearly- 
impenetrable partition with "D" and "N" boundary condi- 
tions enforced at the left and right segment ends, respectively, 
Eq. (|36[l . in units of bulk Fermi energy as a function of 
dimensionless position of the partition na (Fig. la) in the 
^ na ^ TV + 1/2 range for TV = 24 fermions. Bottom: 
corresponding dimensionless Casimir force. 



These observations are illustrated in Fig. 5 where the 
Casimir interaction (top), Eq. (|36[) . is plotted in dimen- 
sionless units of the Fermi energy ef as a function of 
the dimensionless position of the partition na (Fig. la). 
Since the bulk density is (TV + l/2)/L, the variable na 
varies between zero and TV+ 1 /2. Since reflection symme- 
try is lacking, the whole ^ na ^ TV + 1/2 range is shown 
with the same choice of TV = 24 fermions. The bottom 
plot shows the corresponding dimensionless Casimir force 
as a function of na. We notice that except for near a "N" 
wall na = TV + 1/2, the energy and force plots are nearly 
identical to their counterparts, Figs. 2 and 3, with all 
boundaries subject to the "D" wall boundary conditions. 
On the other hand, the immediate vicinity of the right 
boundary, na = TV + 1/2, is similar to the behavior near 
the "N" end, na = 0, in Fig. 4. 

For the cases of the mixed boundary conditions the 
macroscopic limit of the upper envelope of the Casimir 
interaction continues to be given by Eq. ([3TJ]) . Since this 
is the locus of the points of resonant tunneling, then we 
can argue that for a "perfectly transparent" partition and 
any combination of " D" or " N" wall boundary conditions 
enforced at the segment ends, the Casimir interaction is 
given by Ea.pD]). 



IV. DISCUSSION 

In this paper we employed a model of spinless free 
fermions in one dimension to complement the macro- 
scopic approach to the Casimir effect where it is de- 
rived from a regularized low-energy harmonic field the- 
ory. In the segment geometry with like (either DD or 
NN) or periodic boundary conditions when the discrete- 



ness of the underlying particles is irrelevant, we found 
an attractive Casimir interaction of the — nhc/ (2Aa) or 
—nhc/(6a) form, respectively, which is in agreement with 
well-known results of the conformal field theory [22j ■ Ad- 
ditionally we found an attractive Casimir interaction of 
the — irhc/(24a) form for the case of unlike (DN or ND) 
boundary conditions. 

When a nearly-impenetrable partition is inserted in- 
side the segment, a redistribution of the particles takes 
place so as to minimize the system energy; the discrete- 
ness of the fermions plays a major role in determining 
the Casimir force exerted on the partition. However if 
the particle discreteness is ignored, then for any com- 
bination of the D and N boundary conditions enforced 
at the segment ends, the lower envelope of the Casimir 
interaction is found to be universally given by Eq.Q. 

Particle discreteness gives rise to new effects beyond 
the conventional Casimir force, as first noted by Volovik 
[Tl| . However, we disagree with his conclusion that 
for the example he proposed (and we worked out) "the 
Casimir vacuum pressure is not within the responsibil- 
ity of the effective theory, and the microscopic (trans- 
Planckian) physics must be evoked" [26[. This is demon- 
strated by showing how the effects of the particle dis- 
creteness can be understood macroscopically; the result, 
Eq. (|12p , is in agreement with microscopic analysis of 
spinless free fermions. 

The arguments that led to Eq. (fT2"|) highlight the as- 
pects of the Casimir effect in one-dimensional quantum 
liquids which are obscured in the exact microscopical 
treatment of non-interacting fermions: 

First, the effect of particle discreteness accumulated 
in the (na - M - ipf term of Eqs.©, dTUJ) and ||T2J) is 
of nearly-classical origin because the magnitude of the 
effect does not explicitly depend on Planck's constant 
h. This fact is not obvious from the exact analysis of 
free fermions because the latter system does not have a 
classical limit; the sound velocity c in a free-fermion gas 
is determined by Planck's constant H. The only quantity 
which can be due to quantum effects is the phase of the 
Casimir oscillation set by the parameter if. 

Second, the oscillatory part of the Casimir effect is 
not an artifact of the fermionic model, as the low-energy 
theory ^ is applicable both to interacting fermions and 
repulsive bosons. 

Third, the dependence of the effect on the square of 
na — M — <p in Eqs. (f26|) and (f34|) directly demonstrates 
that we are dealing with a phenomenon which is within 
the reach of the harmonic theory. Moreover all the pa- 
rameters entering Eq. (|12|) , even the particle density n 
(which here arises from the reciprocal of the interparticle 
spacing rather than the wavevector cutoff), are macro- 
scopic. Thus there is no need to invoke microscopic 
physics to understand the effect. The only exception to 
this might involve the parameter (p which our macro- 
scopic analysis did not determine. However there is an 
argument that concludes that for arbitrary harmonic liq- 
uid in the DDD geometry one has ip = 1/2 in the L — > oo 
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limit [131 ]. The same value was found rigorously in this 
work for non-interacting fermions and an arbitrary re- 
lationship between a and L. Although more work is 
needed, it seems plausible that for an arbitrary harmonic 
liquid occupying the Dirichlct segment of length L with a 
nearly-impenetrable partition at x = a one has ip = 1/2. 

We see that the continuum contributions to the 
Casimir effect agree with corresponding regularization- 
based conclusions, Eq. ([4]) ; the oscillatory part of the ef- 
fect, up to its phase, is common to all the geometries 
studied, and is correctly described by the macroscopic 
result (PI) . 

To summarize, we have demonstrated that in the 
macroscopic limit, the results of the rigorous analysis of 
the Casimir effect in a one-dimensional gas of spinless 
free fermions are in agreement with the predictions of the 
bosonic low-energy field theory ([3]) if the walls are of the 
Dirichlet type and if the particle discreteness is properly 
accounted for. The latter implements the physical fact of 
near impenetrability which means the number of parti- 
cles captured on either side of the partition is quantized. 
The difference between Eqs.(j4|) and (fT2|) is due to the par- 
ticle discreteness effect; the comparison of Eqs. (fT2")) and 



(|26| (describing specific model under the same physical 
conditions) shows that the Casimir effect for a fermionic 
system is properly described by a macroscopic theory. 

Although it seems unlikely that the Neumann case will 
be any different, it would still be interesting to establish 
what kind of the wavefunction boundary condition cor- 
responds to the Neumann boundary condition imposed 
on the collective variables, and to repeat our analysis for 
this case. This is important because (ignoring the effects 
of the particle discreteness) all the examples we studied 
exhibited only attractive Casimir interaction. Addition- 
ally, it would be very useful to be able to derive (rather 
than guess) macroscopic boundary conditions starting 
from microscopic theory. 
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